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1 Introduction 



The experimental data from DESY HERA on the structure function F2 of deep-inelastic 
scattering (DIS) [T]|2]IMll5]l6]mM and its derivatives (9^2/(9 In |il[6l[T5] 

and 91nF2/91n(l/x) [T5|16fl7|ITS] bring us into a very interesting kinematic range for 
testing theoretical ideas on the behavior of quarks and gluons carrying a very small frac- 
tion of the proton's momentum, the so-called small- a; region. In this limit, one expects 
that the conventional treatment based on the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equations [T9|20|21||22] does not account for contributions to the cross section 
which are leading in a,, ln(l/a;); moreover, the parton density functions (PDFs), in par- 
ticular the one of the gluon, become large, and the need arises to apply a high-density 
formulation of QCD. 

However, reasonable agreement between HERA data and the next-to-leading-order (NLO) 
approximation of perturbative QCD has been observed for ^2 GeV^ (see reviews in 
Ref. [23f2^ and references cited therein) indicating that perturbative QCD can describe 
the evolution of F2 and its derivatives down to very small values, traditionally char- 
acterized by soft processes. 

The standard program to study the x dependence of quark and gluon PDFs is to compare 
the numerical solutions of the DGLAP equations with the data and so to fit the parameters 
of the X profiles of the PDFs at some initial factorization scale Qq and the asymptotic 
scale parameter A. However, for analyzing exclusively the small-x region, there is the 
alternative of doing a simpler analysis by using some of the existing analytical solutions 
of the DGLAP equations in the small-x limit [25|26|27ll28] . This was done in Ref. [25] , 
where it was pointed out that the small-x data from HERA can be interpreted in terms 
of the so-called double-asymptotic-scaling (DAS) phenomenon related to the asymptotic 
behavior of the DGLAP evolution discovered in Ref. [29] many years ago. 

The study of Ref. [25] was extended in Refs. [25|27|28] to include the subasymptotic 
part of the evolution. This led to predictions [271128] of the small-x asymptotic PDF 
forms in the framework of DGLAP dynamics starting at some initial value Qq with flat x 
distributions: 

xfaix, Ql) = Aa (a = q, g), (1) 

where fa{x, Q^) are the PDFs and Aa are unknown constants to be determined from the 
data. We refer to the approach of Refs. [2B|27|28j as generalized DAS approximation. In 
this approach, the fiat initial conditions in Eq. ([T]) play the basic role of the singular parts 
of the anomalous dimensions by determining the small-x asymptotics, as in the standard 
DAS case, while the contributions from the finite parts of the anomalous dimensions and 
from the Wilson coefficients can be considered as subasymptotic corrections, which are, 
however, important for better agreement with the experimental data. In the present 
paper, similarly to Refs. [25|26|I27|I28] . we neglect the contribution from the non-singlet 
quark component. 

The use of the fiat initial condition given in Eq. is supported by the actual experimental 
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situation: small-Q^ data [ilim][T5ll3nll3T]l32] are well described for Q'^ < 0.4 GeV^ by 
Regge theory with Pomeron intercept ap{0) = 1 + Ap = 1.08 (see Ref. [33] and references 
cited therein), close to the standard one, ap{0) = 1. The small rise of the HERA data 
[HI5|11II13|[T5] at small values of can be explained, for instance, by contributions of 
higher- twist operators [28] . 

The purpose of this Letter is to compare the predictions for the structure function 
F2{x,Q'^) and its slope 

91nF2/c)ln(l/a;) from the generalized DAS approach with HI and ZEUS experimen- 
tal data plMEll5llWril8MT0in^ Detailed inspection of the HI data 
points [^]IU1ITB] reveals that, in the ranges x < 0.01 and Q^^2 GeV^, they exhibit a 
power-like behaviour of the form 

F2(x,g2) = Cx"^(«'), (2) 

where the slope A(Q^) is, to good approximation, independent of x and scales logarithmi- 
cally with asA(Q2) = aln(QVA2). A fit yields C ^ 0.18, a ^ 0.048, andA = 292 MeV 
[16j. The linear rise of X{Q'^) with InQ^ is also indicated in Figs. |2]and[3| to be discussed 
below. 

The rise of X{Q^) linearly with InQ^ can be traced to strong nonperturbative physics 
(see Ref. [31] and references cited therein), i.e. A(Q^) ~ l/ag^Q"^). However, the analysis 
of Ref. [22] demonstrated that this rise can be explained naturally in the framework of 
perturbative QCD (see also Section [3]). 

The HI and ZEUS Collaborations [T5|17fl8j also presented new data for A(Q^) at quite 
small values of Q^. As may be seen from Fig. 8 of Ref. [T5|, the ZEUS value for X{Q^) 
is consistent with a constant of about 0.1 at Q^^O.6 GeV^, as is expected under the 
assumption of single-soft-Pomeron exchange within the framework of Regge phenomenol- 
ogy- 
It is interesting to extend the analysis of Ref. [35] to the small-Q^ range with the help 
of the well-known infrared modifications of the strong-coupling constant. We shall adopt 
the "frozen" and analytic ^7] versions. 

This paper is organized as follows. Section[2]contains basic formulae for the structure func- 
tion F2 and its slope 91n F2/51n(l/a;) in the generalized DAS approximation [?r|28|l35] . 
which are needed for the present study. In Section |3| we compare our results on F2 and 
d\nF2/d\n(l/x) with experimental data. Our conclusions may be foimd in Section [ij 



2 Generalized DAS approach 

The flat initial conditions in Eq. ([T| correspond to the case when the PDFs tend to 
constants as a; at some initial value Qq. The main ingredients of the results at the 
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leading order (LO) [27f28] include the following]^ Both, the gluon and quark-singlet PDFs 
are presented in terms of two components ("+" and "— "), 

F^ix, Q2) = exf,{x, Q2), /^(a;, Q^) = Q^) + f-{x, Q^) (a = q,g), (3) 

which are obtained from the analytic Q^-dependent expressions of the corresponding ( "+" 
and "— ") PDF moments. Here, e = Cl2i=i^'i)/f average charge square and / is 

the number of active quark flavors. The small-x asymptotic results for the PDFs are 

xf:{x,Q') = { + ^^.) p/iWe-^+«^ + 0(p), /;(x,Q2) = -^|^/+(x,Q^), 

xf-{x,Q') = V'-^'^^ + f;{x,Q') = -^/-(x,Q2), (4) 

where d+{l) = 1 + 20//(27/?o) and c/_(l) = 16//(27/?o) are the regular parts of the 
anomalous dimensions d^{n) and d_{n), respectively, in the limit n — >■ l|^ Here, n is the 
variable in Mellin space. The functions I^, {v = 0, 1) are related to the modified Bessel 
function and the Bessel function J^, by 

f { \ \ hW), if s > 0; . . 

^^(^) = |rv.(za), if.<0. 

The variables s, a, and p are given by 



a 



where (i+ = — 12//?o, is the strong-coupling constant in the LO approximation, 

and /3o is the first term of its (5 function. 

Contrary to the approach of Refs. [25|26|I27|I28] . various groups were able to fit the available 
data using a hard input at small values of x, of the form x^^, with different values A > 
at small and large values of [33 38 39 4 0|41f42p43f44f45p46f47f48j . At small values, 
there are well-known such results |33j. At large values, this is not very surprising for 
the modern HERA data because they cannot distinguish between the behavior based on a 
steep PDF input at quite large values and the steep form acquired after the dynamical 
evolution from a fiat initial condition at quite small values. 



As has been shown in Refs. [27II28] . the x dependencies of F2 and the PDFs given by the 
Bessel-like forms in the generalized DAS approach can mimic power-law shapes over a 
limited region of x and values: 

Q') ~ x-'^t^^^^'\ xUx, Q') ~ a;-^^^"(^-«^). (7) 



1 The NLO results may be found in Refs. |27l28j . 

^We denote the singular and regular parts of a given quantity k(n) in the limit n — > 1 by k{n) and 
/c(n), respectively. 



4 



In the twist-two LO approximation, the effective slopes have the following forms: 



The corresponding NLO expressions and the higher-twist terms may be found in Refs. PTPS] . 

The effective slopes and A^''^ depend on the magnitudes Aa of the initial PDFs and 
also on the chosen input values of Qq and A. To compare with the experimental data, it 
is necessary to use the exact expressions from Eq. ([s]), but for a qualitative analysis one 
can use some appropriate approximations. At large values of Q^, the "— " components of 
the PDFs are negligible, and the dependencies of the slopes on the PDFs disappear. In 
this case, the asymptotic behaviors of the slopes are given by the following expressions]^ 

where the symbol ~ marks the approximation obtained from the expansion of the usual 
and modified Bessel functions in Eq. (|5]). One can see from Eq. ^ that the gluon effective 
slope X'^J^'^^ is larger than the quark one A^*^'*^*^, which is in excellent agreement with global 
analyses [23ll2l] . 



3 Comparison with experimental data 



Using the results of the previous sections, we analyze HERA data for the structure function 
F2 and its slope 51nF2/c}ln(l/x) at small x values from the HI and ZEUS Collaborations 
[1,2,3,4,5. 6,7. 8. 9|lUflllll2|ll3fl4lll5|ll6fl7lll8j The experimental results for the x dependence 
of F2 in bins of are shown in Fig. [l| while the dependence of A|?(a;,(5^) for an 
average small-x value of 10~^ are shown in Figs. [2]and|3] 

In order to keep the analysis as simple as possible, we fix / = 4 and a^^(M|) = 0.1166, 



so that A^ = 284 MeV and AJ^^^q = 112 MeV, in agreement with the more recent ZEUS 
results [13]. We fit the combined HI and ZEUS data on F2 



mmmmmmmmmmmm 



at LO and NLO imposing two different cuts on Q^, namely > 1.5 GeV^ and > 



0.5 GeV^. The resulting values for Ag, Ag, and Ql are collected in Table [I] together with 
the values of P^r data point (x^/n.d.f.) achieved. In Fig. [T| the HI and ZEUS data 
on F2, which come as x distributions in bins of Q^, are compared with the NLO result 
obtained with the cut > 0.5 GeV^. Furthermore, the dependence of A|^(x,(5^) as 
determined by HI and ZEUS at an average small-x value of 10~^ is confronted with the 
resuh of the NLO fit for > 0.5 GeV^ in Fig. 12} 



■^The asymptotic formulae given in Eq. ^ work quite well at any values > Qf^ because at = Ql 
the values of A^^ and are equal to zero. The use of the approximations in Eq. (pi) instead of the exact 
results given in Eq. ([s]) underestimates (overestimates) the gluon (quark) slope at > Ql only shghtly. 
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Table 1: Results of the LO and NLO fits to HI and ZEUS data 
[T]|21l3P1l5ll5]|7ll5]|9][TUlllllll2lll3lll4] for different small-Q^ cuts. 



UBaaaiaBaaiBiiMiiBiBiH 





A, 




Ql [GeV^] 


xVn.d.f. 


LO 

LO an. 
LO fr. 


0.784 ±0.016 
0.932 ±0.017 
1.022 ±0.018 


0.801 ±0.019 
0.707 ±0.020 
0.650 ±0.020 


0.304 ± 0.003 
0.339 ±0.003 
0.356 ±0.003 


754/609 
632/609 
547/609 


NLO 
NLO an. 
NLO fr. 


-0.200 ±0.011 
0.310 ±0.013 
0.180 ±0.012 


0.903 ±0.021 
0.640 ±0.022 
0.780 ±0.022 


0.495 ±0.006 
0.702 ±0.008 
0.661 ±0.007 


798/609 
655/609 
669/609 


g2 > 0.5 GeV2 
LO 

LO an. 
LO fr. 


0.641 ±0.010 
0.846 ±0.010 
1.127± 0.011 


0.937 ±0.012 
0.771 ±0.013 
0.534 ±0.015 


0.295 ±0.003 
0.328 ±0.003 
0.358 ±0.003 


1090/662 
803/662 
679/662 


NLO 
NLO an. 
NLO fr. 


-0.192 ±0.006 
0.281 ±0.008 
0.205 ±0.007 


1.087 ±0.012 
0.634 ±0.016 
0.650 ±0.016 


0.478 ±0.006 
0.680 ±0.007 
0.589 ±0.006 


1229/662 
633/662 
670/662 



Because the twist-two approximation is only reasonable at ^2.5 GeV^ [28], as may be 
seen from Fig. [l| some theoretical improvements are necessary for smaller values. In 
Ref . [28] , the higher- twist corrections through twist six were added to find good agreement 
for (5^^0.5 GeV^. However, the twist-four and twist-six terms increase the number of 
parameters, which become strongly correlated. 

Here, we investigate an alternative possibility, namely to modify the strong-coupling con- 
stant in the infrared region. Specifically, we consider two modifications, which effectively 
increase the argument of the strong-coupling constant at small values, in accordance 
with Refs. |49f50f51|l52f53f54j . In the first case, which is more phenomenological, we intro- 
duce a freezing of the strong-coupling constant by changing its argument as (5^±M^, 
where Mp is the rho- meson mass [M] • Thus, in the formulae of Section [2] and their NLO 
generalizations [271)28] . we introduce the following replacement 

- "L(Q') = "UQ' + K) (i = LO,MS), (10) 

where a^*~'((5^) and q;^^((5^) have the canonical forms dictated by the renormalization 
group. 

The second possibihty is based on the idea by Shirkov and Solovtsov [37f55] (see also the 
recent reviews in Refs. [56|57f58j and the references cited therein) regarding the analyticity 
of the strong-coupling constant that leads to an additional power dependence. In this 
case, the one-loop and two-loop coupling constants a^°(Q^) and a^^(Q^) appearing in 
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the formulae of the previous sections and their NLO generahzations are to be replaced as 

1 A2 

a, [Q j ftan ) = [Q 



«r (Q^) - c^Tm = af'm - ^ + ..., (11) 

where the ellipsis stands for cut terms which give negligible contributions. 

We thus repeat the LO and NLO fits discussed above using in turn the frozen and analytic 



versions of the strong-coupling constant according to the replacements of Eqs. (10) and 
(11), respectively. The results for F2 are included in Table [l] and Fig. [l] and those for 
A^(x, Q^) in Fig. 2j Figure [l| as also Fig. 5 in Ref. |28], only covers the small-Q^ region, 
Q'^ < 9.22 GeV^, because the three NLO predictions are hardly distinguishable at larger 
values of Q^. 

From Fig. [T], we observe that the fits based on the frozen and analytic strong-coupling 
constants are very similar and describe the data in the small-Q^ range significantly better 
than the canonical fit. This is also reflected in the values of x^/n.d.f. listed in Table [!} 
The improvement is especially striking at NLO if data with very small values, with 
> 0.5 GeV^, are included in the fits. Then x^/n.d.f. is almost reduced by a factor of two 
to assume values close to unity when the canonical version of the strong-coupling constant 
is replaced by the frozen or analytic ones. The situation is very similar to the case when 
the higher-twist corrections according to the renormalon model are included f2H|. In order 
to illustrate this, we display in Fig. [TJalso the results obtained at NLO in the renormalon 
model of higher-twist terms, which are taken from Fig. 5 in Ref. [28]. We see that the 
latter describe the experimental data shghly better for 0.65 GeV^ ^2.0 GeV^ than 
the results obtained here, which is also reflected in the values of x^/n.d.f. achieved, namely 
X^/n.d.f. = 565/658 = 0.86 for renormalon improvement versus x^/n-d.f. = 633/662 = 
0.96 and 670/662 = 1.01 for analytic and frozen strong-coupling constants, respectively. 
However, one should bear in mind that this improvement happens at the expense of 
introducing four additional adjustable parameters. 

Figure [2] nicely demonstrates that the theoretical description of the small-Q^ ZEUS data 
on A^(a;, Q"^) by NLO QCD is significantly improved by implementing the frozen and ana- 
lytic strong-coupling constants. Again, these two alternatives lead to very similar results. 
For comparison, the linear rise of A|^(x, Q^) with In as described by Eq. ^ is also indi- 
cated in Fig. |2} For comparison, we display in Fig. |2] also the results obtained by Kaidalov 
et al. |17] and by Donnachie and Landshoff [H] adopting phenomenological models based 
on Regge theory. While they yield an improved description of the experimental data for 
^ 0.4 GeV^, the agreement generally worsens in the range 2 GeV^ ^ ^ 8 GeV^. 



As may be seen from Table 



(see also Ref. 



Figure 



|2pi 



the three NLO fits for F2{x,Q^) yield Ql ^ 0.5-0.7 GeV^ 
shows that the conventional NLO fit yields A^(x, Ql) = as 



suggested by Eq. (T|. The replacements of Eqs. ([To ) and ( 11 ) raise the value of A^(x, Q 
In fact, the results for X'p {x,Q'^) obtained with t " 



le frozen and analytic versions of the 
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strong-coupling constant agree much better with the ZEUS data at GeV^. Nev- 

ertheless, for < 0.5 GeV^, there is still some disagreement with the data, which needs 
additional investigation. 

In Fig. |2| the NLO results for Xf^{x, Q^) are evaluated at x = 10 ^. In Fig. |3| we study 
the variation with x in the range 10~^ < x < 10~^. For simplicity, we only do this for the 
case of the frozen strong-coupling constant; the result for the analytic one would be very 
similar. We observe good agreement between the experimental data and the generalized 
DAS approach for a broad range of small-x values. At small values, A|^(x,(5^) is 
practically independent of x, which is because the variable p defined in Eq. ^ takes rather 
small values there. At large values, the x dependence of A^(x,(5^) is rather strong. 
However, it is well known that the boundaries and mean values of the experimental x 
ranges [16] increase proportionally with Q"^, which is related to the kinematical restrictions 
in the HERA experiments, namely x ~ 10~^ x (see Refs. [If6|13fl4] and, for example. 
Fig. 1 of Ref. [I5]). From Fig. |3| we see that the HERA data are close to A|f(x,(5^) 
at X ~ 10-^-10-5 for = 4 GeV^ and at x ~ lO'^ for = iqo GeV^. Indeed, the 
correlations between x and of the form Xeg = a x 10~^ x with a = 0.1 and 1 lead 
to a modification of the evolution which starts to resemble InQ^, rather than InlnQ^ 
as is standard 135). 



4 Conclusions 



We studied the dependence of the structure function F2 and the slope = d\nF2/d\n{l/x) 
at small x values in the framework of perturbative QCD. Our twist-two results are 
in very good agreement with HERA data [TllMllllSllGllTllSlI^^ at 
Q^^2.5 GeV^, where perturbation theory is applicable. The applications of the frozen 
and analytic versions of the strong-coupling constants, a^^{Q'^) and a^{Q^), significantly 
improve the agreement with the HERA data p!P|5|lf5|5|7|5|5|TU|ll|12|ll3|14|15lll6|17|18j 
for both the structure function F2 and its slope A^(x, Q^) for small values, ^ 0.5 GeV^. 
The results obtained with these infrared-modified strong- coupling constants and also those 
based on the renormalon model with higher-twist terms incorporated, which were consid- 
ered in Ref. [28j, are very similar numerically. 



As a next step of our investigations, we plan to fit the HERA data [T|2|3| 
for F2{x,Q'^) using alternative analytic versions of the strong-coupling constant 



see. 



for 



example, the recent reviews in Refs. [^571158] ). to find out if the theoretical description 
of the slope A^ can be further improved at small values. 
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Figure 1: x dependence of F2{x, Q^) in bins of Q^. The experimental data from HI (open 
points) and ZEUS (solid points) are compared with the NLO fits for > 0.5 GeV^ im- 
plemented with the canonical (solid lines), frozen (dot-dashed lines), and analytic (dashed 
lines) versions of the strong-coupling constant. For comparison, also the results obtained 
in Ref. [28] through a fit based on the renormalon model of higher-twist terms are shown 
(dotted lines). 
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Figure 2: Q'^ dependence of A^(x,Q^) for an average small-x value of x = 10"^. The 
experimental data from HI (open points) and ZEUS (solid points) are compared with the 
NLO fits for > 0.5 GeV^ implemented with the canonical (solid line), frozen (dot- 
dashed line), and analytic (dashed line) versions of the strong-coupling constant. The 
linear rise of A|?(x,Q^) with InQ^ as described by Eq. ([2) is indicated by the straight 
dashed line. For comparison, also the results obtained in the phenomenological models 
by Kaidalov et al. [57] (dash-dash-dotted line) and by Donnachie and Landshoff 
(dot-dot-dashed line) are shown. 
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Figure 3: Same as Fig. [2| but for the NLO fit implemented with the frozen version of the 
strong-couphng constant and for x in the range 10^^ < x < 10^^. 
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